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ABSTRACT 

We performed 3D numerical simulations of the solar surface wave field for 
the quiet Sun and for three models with different localized sound-speed varia- 
4^ tions in the interior with: (i) deep, (ii) shallow, and (iii) two-layer structures. 

We used simulated data generated by two different codes which use the same 
standard solar model as a background model, but utilize two different integra- 
tion techniques and use different models of stochastic wave excitation. Acoustic 
travel times were measured from all data sets using the time-distance helioseis- 
mology technique and compared with the ray theory predictions, frequently used 
for helioseismic travel-time inversions. It is found that the measured travel-time 
shifts agree well with the ray theory in both cases with and without phase- 
speed filtering for the shallow and deep perturbations. This testing verifies the 
whole measuring-filtering-inversion procedure for sound-speed anomalies inside 
the Sun. It is shown, that the phase-speed filtering, frequently used to im- 
prove the signal-to-noise ratio does not introduce significant systematic errors. 
Results of the sound-speed inversion procedure show good agreement with the 
background sound-speed profiles in all cases. Due to its smoothing nature, the 
inversion procedure overestimates sound speed variations in areas with sharp 
gradients of the sound-speed profile. 

Subject headings: Oscillations, Solar; Magnetohydrodynamics; Sunspots, Mag- 
netic Fields 



1. Introduction 



Both ti me-distance helioseismol ogy (jDuvall et al.lll993l ) and phase- sensitive helioseismic 
holography (IBraun fc Lindseyl 120001 ) measure acoustic travel-times or acoustic wave phase 
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shifts to infer solar interior properties such as sound-speed perturbations and flow fields. 
For time-distance helioseismology inversions, the measured travel-time anomalies associated 
with sunspots and active regions are often i nterpr e ted as caused by solar subsurface sound- 



speed perturbations (e.g. iKosovichev et al.l (120001 ); ICouvidat et al.l (120041 ) ). The travel-time 



inversions for the sound-speed p erturbations are usually don e by using a variational formula 
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20041 ) and provide similar results for sunspot data 



It is important to verify and calibrate the standard helioseismic inversion procedure, 
which can be done by analyzing artificial data obtained from numerical simulations. We ap- 
ply the time-distance measurement and helioseismic inversion procedures to numerically sim- 
ulated stochastic velocity wave fields on the Sun calculated using di fferent simulation codes. 
Both c odes solve the linearized Euler equations. The Cartesian code (jParchevsky fc Kosovichev 
2007al ) uses a finite-difference scheme with impulsive localize d wave sources randomly d is- 
tributed at same depth and random in time. Another code (IHartlep fc Mansourl 120051 ) is 
spherical and utilizes a combination of a spectral method in angular coordinates and B- 
splines in radial direction. Wave sources used in this code are distributed and defined in the 
Fourier space as a combination of spherical harmonics with random amplitude and phase 
and a mean acoustic power spectrum similar to the solar one. 

Our goal is to examine how well the time-distance measurements and sound-speed in- 
version results based on acoustic ray-path approximation agree with the travel-time shifts 
expected from the acoustic ray theory and sound-speed variations of the background model. 
Similar testing for the travel times measured by the acoustic holography method and sound - 
speed inversions based on the Born approximation has been published by iBirch et al.l ( 120111 ) . 
We introduce the models and numerical simulations in §2. In §3 we describe details of the 
time-distance measurement techniques with and without phase-speed filtering. In §4 we 
present results of sound-speed inversions for different profiles of the background sound speed 
variations using the artificial data, obtained by different codes. Conclusions follow in §5. 



2. Artificial data 



2.1. Simulation of local wavefield 



We use the linearized c ode for simulations of wave propaga tion in the Cartesian geometry 
in the Sun, described by iParchevsky fc Kosovichevl (J2007aj). A finite difference scheme, 
adapted to the non-uniform vertical grid, is used for the spatial discretization. A strong- 
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stability-preserving Runge-Kutta algorithm is used for time advancing ( jShul 120021 ). We 
simulated the wave field generated by localized sources with amplitude A, added in the right 
hand side of the momentum equation as a vertical force 



S z (r,t)=AG(r)F(t), 



(1) 



randomly distributed in a shallow layer 100 km below the photosphere. The spatial G(r) 
and temporal F(t) behavior of the wave source is modeled by functions 



G(r) 



1 - 



H 2 
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where R src and H src are the source radius and height respectively, x src , y src , 
coordinates of the source origin. Parameter r is given by equation 



and 



u (t - 1 ) 
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where uq is the central source frequency, to is the moment of the source initiation. In 
our simulations we chose H src = 0.1 Mm, R src = 0.3 Mm. This source model provides 
the wave spectrum, which closely resembles the solar spectrum. It has a peak near the 
central frequency cup an d spreads over a broad frequency interval. The source spectrum is 
dParchevskv et al.ll2Q08h : 



F(t)e 



-•wot 



dt 



UJ 



e o . 



(5) 



The acoustic wave field, generated by this code with described distribution of wave sources, 
has a power spectrum that is very similar to the spectrum obtained from SOHO/MDI high- 
resolution observations (Fig. [1]). 

Great efforts were made to develop the top boundary condition which mimics non- 
reflective properties of the solar atmosphere. Waves with frequencies below the acoustic cut- 
off frequency a re reflected fr om the sub-photospheric layer with strong density and pressure 
gradients (e.g. iGoughl (119931 )). The waves with frequencies higher than the acoustic cut-off 
frequency, pass through and propagate to the chromosphere and escape to the corona. For 
simulation of this property, we set a wave absorbing layer, based on Perfectly Matched Layer 
(PML) technique, in low chromosphere, about 500 km above the photosphere. The waves 
with the frequencies higher than the acoustic cut-off frequency are absorbed by this layer 
and do not pollute the computational domain. The higher layers are not considered because 



-4- 



they are not essential for the current helioseismology study. A similar absorbing layer was 
set at the bottom of the domain. The lateral boundary conditions were chosen periodic. 

We built a background plai n-parallel model of the Sun by smoothly joining the model 
of the chromosphere ("VALQ of IVernazza et al.l (119811 ) with the standard solar model S of 
Christensen-Dalsgaard et al.l (119961 ). To eliminate the convective instability in the subsurface 
superadiabatic region, the background model was stabilized aga inst convection by setting 
the neg ative entropy gradient in the entropy equation to zero (see iParchevsky fc Kosovichev 
j2007al i). 



2.2. Simulations of whole-Sun oscillations 



Si mulations of the global full- Sun acoustic wave field are preformed with a code devel- 
oped by lHartlep fc Mansourl (120051 ). This code solves the linearized propagation of helioseis- 
mic waves throughout the entire solar interior, and has been used in various studies to sim- 
ulate the effects of localized sound speed perturbations, e.g., for testing helioseismic far-side 



imaging by simulati ng the effects of model sunspots on the acoustic field (IHartlep et al.ll2008 



Honidis et al.l 120091). for validat ing time-distance helioseismic measurements of tachocline 
perturb ations (jZhao et al.l 120091 ) , and for studying the effect of localized subsurface pertur- 



bations (IHartlep et al 



20111 ). It models solar acoustic oscillations in a spherical domain by 



using the Euler equations linearized around a stationary, non-rotating, background state, 
and uses several approximations. In particular, perturbations of the gravitational potential 
have been neglected, and the adiabatic approximation has been used. Also, the entropy 
gradient of the background model has been neglected in order to make the linearized equa- 
tions convectively stable. The calculations show that this assumption does not significantly 
change the propagation properties of acoustic waves including their frequencies, except for 
the acoustic cut-off frequency, which is slightly reduced. Non-reflecting boundary conditions 
are used at the upper boundary by means of an absorbing buffer layer with a damping 
coefficient that is zero in the interior and increases smoothly into the buffer layer. 

A random function was added in the mass continuity equation to mimic the excitation 
of acoustic sources, which is a non-linear process not included the in linearized equations. 
This function is non-zero only close below the solar surface. It is random in latitude and 
longitude, and in time, but with a prescribed temporal spectrum approximating the solar 
acoustic power spectrum. 

The equations are solved using a pseudo-spectral metho d with a spherical harmonic dis- 
cretization in the angular directions, and 4th order B-splines ( Loulou et al.|[l~997 ; Kravchenko et al 
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19991 ) for the radial direction. 2/3-dealiasing is used. The radial resolution of the B-spline 



method is varied proportionally to the local speed of sound so that the method's resolution 
in greater near the surface (where the sound speed is small), and coarser in the deep interior 
(where the sound speed is large). The simulation in this study resolves spherical harmonics of 
angular de gree between to 170, and uses 350 B-splines in the radial direction. A staggered 



Yed (119661 ) scheme is used for time integration, with a time step of 1 seconds. 



There are two main differences between the Cartesian (CC) and the spherical (SC) 
codes: (i) CC uses finite difference spatial scheme, SC uses decomposition of the wave field on 
spherical harmonics in the horizontal direction, (ii) CC uses vertical dipole sources localized 
in space and time and defined in the configurational space, SC code uses sources defined 
in the Fourier space as a composition of spherical harmonics with random parameters. So, 
in the configurational space, sources are not localized, but rather smeared over the whole 
solar surface. These two different approaches to the simulation of the wavefield in the Sun 
produce the acoustic wave fields with similar parameters. 



3. Travel-time measurement technique 

Travel-time measurements are based on fitting of the cross-covariance function 

C(n, r 2 , t) = ^f 4>(ri, t')4>(r 2 , t' + t)t' (6) 

between two source and destination points with horizontal radius- vectors ri and r 2 , where 
T is the total time of observation, t is the time lag between two signals, and is the vertical 
(z-component) of the velocity perturbation in the wave. An /-mode filter is applied to 
4>{ri,t) and 0(r 2 ,t) to filter out the /-mode. Then, to increase the signal-to-noise ratio, the 
point-to-point cross-covariance functions are averaged over the annulus with internal and 
external radii A min , A max , centered at the source. To do this, the annulus is divided along 
the radius into pixel-wide sub-annuli and point-to-annulus cross-covariance functions are 
calculated for each sub-annulus by averaging the point-to-point cross-covariance functions 
along each sub-annulus. Then, the individual point-to-annulus cross-covariance functions 
are shifted in time in such a way, that they are centered on the specific reference time, 
which corresponds to the travel time between the source and the sub-annulus with radius 
l/2(A min + A max ). To produce the final point-to-annulus cross-covariance function, the time 
shifted cross-covariance functions for individual sub-annuli are averaged. 
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Fig. 1. — Oscillation power spectra from high-resolution SOHO/MDI observations (top), and 
from the two numerical simulations of the quiet Sun in the Cartesian geometry (middle), 
and in the spherical geometry (bottom) as functions of the phase speed, u/k, and cyclic 
frequency, v = u/2n. The leftmost ridge corresponds to the /-mode (surface gravity wave), 
and other ridges represent acoustic (p) modes. The vertical lines indicate where the phase- 
speed filter allows the passage of a half of the power, and the horizontal lines show the 
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Kosovichev Duvalll (119971 ) proposed to fit the cross- covariance function by the Gabor's 
wavelet 

G(A,uj ,5uj,T p ,T g ;t) = Acos(u (t - r p )) exp ( — (t-T g ) 2 j , (7) 



where A is the amplitude, uq is the central frequency, 5co is the frequency width, t p is the 
phase travel-time, and r g is the group travel-time. Equation fl7J) can be considered as a 
definit ion of the phase- and group travel times. This definition is not unique. iGizon fc Birch 
( 120021 ) proposed a definition of the wave packet travel time based on geophysical studies. 
Later, IGizon fc Birch! ( 12004 ) proposed another definition of the travel-t ime linear in the 
cross- covariance. This approach was described by iBirch et al.l (120091 1201 ll ) 



Usually, before the averaging, 4>( r, t) is filtered in t he Fourier domain using a phase- 
speed filter with the Gaussian profile (IDuvall et al.lll997l ). Such filtering selects only waves 
which travel with phase speed in a given range, or, in other words, waves which travel-time 
between the source and destination lies in some range which depends on the width of the 
pase-speed filter. This further improves the signal-to-noise ratio. The central phase-speed 
and the width of the phase-speed filter depend on the internal and external radii A min , A max 
of the annulus and must be carefully tuned out. Here we perform the time-distance analysis 
of the simulation data with and without the phase-speed filtering procedure in order to test 
the accuracy of the filtering. 



4. Results and discussion 

In this section we present results of numerical simulations, time-distance analysis of the 
artificial data, and results of the helioseismic inversion of the sound speed profiles for differ- 
ent models of the sub-photospheric sound speed perturbations. To simulate local subsurface 
sound-speed structures, two types of axially symmetric perturbations of the background 
sound speed were introduced: constant-sign perturbation (two models with different param- 
eters) and two- layer perturbation (one model), where the background sound speed pertur- 
bation changes sign at some depth. The spatial distribution of constant-sign sound-speed 
perturbations is given by equation 



5c 2 (r,z)/c 2 (z) = eexp(-r 2 /R 2 spt )exp(-z 2 /D 



2 > 
sptJ 



where r is the horizontal distance from the symmetry axis of the perturbation, z is the 
depth relative to the photosphere, c (z) is the sound speed in the quiet Sun model, R spt is 
the horizontal scale of the sound speed perturbation, D spt is the depth scale, and e is the 
fractional amplitude of the perturbation in the sound spe ed frame. We adop ted two different 



sets of parameters for the constant-sign perturbation from lBirch et al.l (120091 ) : "shallow" case 
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Fig. 2. — (a) Example of the cross-covariance function (time-distance diagram) from the 
numerical simulations, using the Cartesian geometry, averaged over the selected central area. 
White dashed curve shows the Gabor's wavelet fitted phase-travel times, (b) Comparison 
of the phase travel times measured from the three simulation models: quiet Sun (thin solid 
curve), cases with shallow (dashed curve) and deep (solid curve) perturbations, (c) Travel 
time shifts obtained without phase-speed filtering (triangles) and with phase-speed filtering 
(squares) from the simulation data relative to the quiet Sun. Dashed and solid curves 
represent time shifts calculated from the acoustic ray theory for the shallow and deep models 
respectively. 




Fig. 3. — Vertical (top panel) and horizontal (bottom panel) slices of the averaging kernels 
used for the helioseismic inversions of simulated data for the shallow, deep and two-layer 
models, described in the text. The solid, dotted, dashed, dash-dotted, and dash-dot-dotted 
curves correspond to the depth ranges of 0-1, 1-2, 2-4, 4-6, 6-8 Mm respectively. 




Fig. 4. — Model and reconstructed sound speed profiles for the shallow model. The top left 
panel shows vertical cut of the sound speed perturbation, used as the background state in 
the shallow model. The top right panel shows the sound-speed perturbation obtained by 
inversion of travel times. The bottom left panel shows the sound-speed perturbation at the 
axis of the shallow model as a function of depth. The dash dotted curve represents the model 
profile given by equation (JSJ), the solid step curve shows the mean reconstructed profile. The 
error bars represent the RMS deviations of the reconstructed sound speed perturbation in 
horizontal plane. The bottom right panel represents the ID horizontal profiles of the sound 
speed perturbation averaged over depth ranges 0-0.5, 1.5-3, 5-8, and 8-11 Mm, shown by 
solid, dashed, dot-dashed, and dotted curves respectively. The thick curves show profiles in 
the background model, the thin curves represent inversion results. 
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(a) with D = 1 Mm, R = 20 Mm, and e = 0.15, and a "deep" case (b) with D — 10 Mm, 
R = 20 Mm, and e = 0.1. Sound-speed perturbations in both "deep" and "shallow" cases 
are positive. The two-layer sound-speed perturbation is negative in a shallow region below 
the photosphere and positive for depths greater than 3 Mm. The vertical profile of the 
sound speed in the two-layer model was chosen to re semble the initial helioseismic sound- 



speed inversion of a sunspot (jKosovichev et al.ll2000i ). The horizontal profiles in all three 



cases were chosen Gaussian with R = 20 Mm. In addition to these simulations with three 
different models of the background sound-speed perturbation, we performed simulations 
with the horizontally uniform "quiet" solar model. Locations, amplitudes, frequencies, and 
moments of initiation of the wave sources in simulations with the quiet solar model and all 
three perturbed models were the same. 

We performed 3D numerical simulations in the Cartesian geometry for all three models 
of sound-speed perturbations: (i) "deep", (ii) "shallow", and (iii) two-layer perturbation. 
Simulations were carried out in a box of size 96 x 96 x 36.8 Mm 3 using grid size of 640 x 640 x 92. 
The horizontal grid is uniform (Ax = Ay = 150 km), and the vertical grid is non- uniform 
with spacing Az = 50 km above the photosphere and gradually increasing to Az = 1.76 Mm 
at the bottom of the of the computational domain. The time step was chosen to be 0.5 s 
to comply the CFL stability condition of the explicit numerical scheme. In all three models 
localized stochastic sources of the vertical component of force with random amplitudes and 
frequencies were uniformly distributed at random locations at the depth of 100 km below 
the photosphere. We use 512-min time sequences with a 1-min cadence of the photospheric 
vertical velocity in a 15 x 15 Mm 2 area for the time-distance analysis. 

The oscillation power spectra obtained from the MDI and the simulations data are 
plotted in Figured] as functions of the oscillation cyclic frequency, v = co/2ti, and the phase 
speed, u/k, where k is the wavenumber, u is the angular frequency. Position of the /- and 
p-mode ridges in all three spectra are nearly identical. The ridge width in the MDI power 
spectrum and power spectrum obtained from simulations in the Cartesian geometry is the 
same. The power distribution of the simulated wavefield in the Cartesian geometry is a bit 
stronger in the high-frequency and low phase-speed region, and a bit weaker in the low- 
frequency and high-phase-speed region. The ridge width of the power spectrum obtained 
from the simulations in the spherical geometry is smaller than in the MDI spectrum (probably 
due to lower dissipation in the sponge region). The simulated wavefield obtained from the 
spherical simulations exhibits excess of power at frequencies below 2.5 mHz and lack of power 
at frequencies above 4 mHz due to the reduced acoustic cut-off frequency caused by the non- 
inclusion of the entropy gradient. However, from our experiments of applying frequency 
filters, we find that these differences between the observed and simulated power spectra do 
not significantly alter the travel-time measurements in this work. The differences in power 




Fig. 5. — Model and reconstructed sound speed profiles for the deep model. Position of 
panels and curve notations are the same as in Figure 0] 
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spectra from simulations in the Cartesian and spherical geometries are due to the different 
formulations of wave source excitation used in these simulations. 

The cross- covariance function (time-distance diagram), obtained from the numerical 
simulations with the quiet-Sun background model is shown in panel (a) of Figure [2j The 
cr oss-covariance funct ion was calculated without phase-speed filtering which was introduced 
by lDuvall et al.l ( 119971 ) to improve the signal-to- noise ratio for the function. The fitted phase 
travel times for three simulations: quiet-Sun, and shallow and deep perturbation background 
models are shown in panel (b). A comparison of the phase travel times clearly shows that for 
the cases with both shallow and deep perturbations, the measured travel times are shorter 
than the travel times for the quiet Sun model, calculated from the acoustic ray theory, as 
expected for the positive background sound-speed perturbations. The phase travel time shifts 
measured with and without phase-speed filtering are shown in panel (c). The horizontal 
error bars indicate the annulus radius ranges used in this analy sis. The vertical bar s are 
formal fitting errors. The filtering parameters are adopted from (jCouvidat et al.l 120061 ) and 
illustrated by vertical lines in Figure [TJ Because of the limited size of the simulation box 
only the first few filt e rs are applied in our study. Essentially the same filter parameters were 
used by iBirch et al.l ( 120091 ). Clearly, there is a good agreement between the measurements 
with and without the phase-speed filtering. At most points, the theoretical values fall within 
the scale of error bars of the measured travel times. This analysis provides an important 
validation for the basic technique of time-distance helioseismology. Of course, the acoustic 
ray theory has well-known limitations. It is expected that the theoretical description of 
acoustic travel times can be im proved by accounting for the finite wavelength effec ts, e.g. by 
using the Born approximation (IBirch fc Kosovichevil2000l ; IBirch et al.l 12004 . 120111 ). 



The sound-speed perturbations for the shallow model are presented in Figure HI The 
top panels represent the background model sound-speed perturbation (left) and sound-speed 
perturbation reconstructed from the travel times by the helioseismic inversion (right). Bot- 
tom left and right panels represent the vertical and horizontal profiles of the sound-speed 
perturbation respectively. The reconstructed and model profiles show very good agreement 
above 2 Mm. The inversion results overestimate the sound-speed variations in the depth 
range from 1.5 Mm to 3.0 Mm, because the boundary of the sound-speed perturbation is 
located near the 1.5 Mm. The reconstructed sound-speed profile is smoothed, which causes 
a slight increase of the perturbation depth, which one can expect from the inversion proce- 
dure. Below 3.0 Mm both model and reconstructed sound speed perturbation profiles are 
zero. Similar results for the deep perturbation model are shown in Figure [5J As for the 
shallow case, the reconstructed profiles of the sound-speed perturbation show good agree- 
ment with the sound-speed perturbation profiles in the background model except the depth 
region between 8 Mm and 11 Mm, which is close to the bottom boundary of the sound speed 
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Depth, Mm Distance, Mm 

Fig. 6. — Model and reconstructed sound speed profiles for the two-layer model. Position of 
panels and curve notations are the same as in Figure HI 
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perturbation. 

Results of reconstruction of the sound-speed perturbation in two-layer model are shown 
in Figure O Description of panels is the same as in Figure HJ The left and right top panels 
represent maps of the model and reconstructed sound speed perturbations respectively, the 
bottom left and right panels show vertical and horizontal profiles of the sound-speed per- 
turbatio n. For comparison , we present inversion results obtained from the same simulation 



data by IBirch et al.l (j201l[ ). They are shown by the dashed curve in the left bottom panel. 
Dotted curves show the error range. The reconstructed profiles are in good agreement with 
the model profiles. The depth of the sign inversion of the sound speed perturbation is located 
at the correct depth. Reconstructed amplitudes of the perturbation are close to the model 
ones, though the inversion slightly overestimates the perturbation amplitude in dept range 
from 1.5 Mm to 3.0 Mm, and are noisy in the shallow subphotospheric region (depth range 
from Mm to 0.5 Mm). 

We also performed numerical simulation of the wave field in the two-la yer model of the 



sound speed perturbation in the spherical geometry using the spectral code of lHartlep fc Mansour 



( 120051 ). It is important to test the measurement and inversion procedures on different sets 
of artificial data, generated by the different codes. The spherical code allows to perform the 
testing for deeper perturbations than the Cartesian code, the depth of the perturbation was 
increased up to 20 Mm, when the depth of the original perturbation is 10 Mm. 

Results of the sound speed inversion in the two-layer model from simulations in spherical 
geometry with the spectral code are presented in Figure [71 We see good agreement between 
the background and reconstructed profiles of the sound speed perturbation, though the 
transition between negative and positive perturbations in the inversion results is not as 
steep as in the background model due to smoothing properties of the inversion procedure. 



5. Conclusions 

By using numerical simulations of stochastically excited solar oscillations we have tested 
the measurement and inversion procedures of time-distance helioseismology, based on the 
Gabor wavelet fitting to the cross-covariance function and ray-path sensitivity kernels. We 
have measured the travel time shifts caused by the subsurface sound-speed perturbations 
with and without phase-speed filtering for three different background distribution of the 
sound speed perturbation. We used two sets of simulated data generated by two different 
codes. One code uses finite difference spatial scheme in the Cartesian geometry and dipole 
wave sources localized in space and time. Another code uses spherical harmonic expansion 
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Depth, Mm Distance, Mm 

Fig. 7. — Reconstruction of the sound speed perturbation for the two-layer model calculated 
from simulations in spherical geometry with the spectral code. Panel description is the same 
as in Figure HJ The sound speed perturbation used as a background model in spherical 
simulations is shown by the dash-dotted curve in bottom left panel. The solid, dashed, 
dot-dashed, and dotted curve in the right bottom panel represent horizontal profiles of the 
sound-speed perturbations averaged over depth ranges 0-1, 3-5, 7-10, 13-17 Mm respectively. 
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in spherical geometry and distributed sources defined in the Fourier space. Simulations have 
been done for three types of the background sound speed perturbations: (i) shallow, (ii) 
deep, and (hi) two-layer models. 

We applied the time-distance analysis to both types of simulations and then performed 
the helioseismic inversion of the obtained travel times to reconstruct the sound-speed per- 
turbation profiles. The time-distance analysis of simulation results for deep and shallow 
models show good agreement with the theoretical time shifts calculated from the acoustic 
ray theory for both cases with and without phase speed filtering. The travel time variations 
in the perturbed regions are negative everywhe re and do not exhibit sign inversion at short 



distances as was reported in iBirch et all (120091 ) for the acoustic holography measurements 



The helioseismic inversion of the travel times for the sound-speed show good agreement 
with the background model profiles for all three cases and for both types of artificial data, 
obtained with the Cartesian and spherical codes. The inversion procedure produces smoothed 
results as one can expect, resulting in overestimation of the sound speed perturbation in the 
shallow case near relatively sharp bottom boundary of the background perturbation. For 
the two-layer model the smoothing nature of the inversion procedure results in the less steep 
gradient of the sound speed perturbation near the depth where the perturbation changes 
sign. These systematic effects should be taken into account in similar analysis of the solar 
oscillation data. 

We summarize our conclusions as follows. 



• The measurements of acoustic travel times from numerical simulations of solar os- 
cillation data using the time-distance helioseismology technique clearly demonstrate 
that the travel time shifts caused by subsurface sound-speed perturbations can be well 
inverted by using the acoustic ray approximation. 

• The test measurements with and without the phase-speed filtering procedure demon- 
strate that this procedure used to improve the signal-to-noise ratio does not introduce 
significant systematic errors. 

• Time-distance analysis and inversions of travel times, obtained from simulations per- 
formed by different codes show that the results are not simulation dependant. The 
sound-speed inversion results show good agreement with the background sound-speed 
profiles in all cases. 
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